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Does sub-cluster merging accelerate mass segregation in 
local clusters? 
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ABSTRACT 

The nearest site of massive star formation in Orion is dominated by the Trapezium 
subsystem, with its four OB stars and numerous companions. The question of how 
these stars came to be in such close proximity has implications for our understand- 
ing of massive star formation and early cluster evolution. A p romising route toward 
rapid mass segregation was proposed by 
merger product of faster-evolving sub c 



McMillan et al.| ( |2007| ), who showed that the 
usters can inherit their apparent dynamical 
age from their progenitors. In this paper we briefly consider this process at a size and 
time scale more suited for local and perhaps more typical star formation, with stellar 
numbers from the hundreds to thousands. We find that for reasonable ages and cluster 
sizes, the merger of subclusters can indeed lead to compact configurations of the most 
massive stars, a signal seen both in Nature and in large-scale hydrodynamic simula- 
tions of star formation from collapsing molecular clouds, and that sub-virial initial 
conditions can make an un- merged cluster display a similar type of mass segregation. 
Additionally, we discuss a variation of the mi nimum spanning tree mass-segregation 
technique introduced by Allison et al. (2009a I. 
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1 INTRODUCTION 

The question of how dense, massive groupings like the 
Trapezium are formed has yet to find a convincing answer. 
Ignori ng for the moment the striking multiplicity of the sys- 



(as opposed to gas-driven segregation) is the young age of 
the cluster and the segregation timescale, which taken to- 
gether at face value would suggest that the cluster is not 



terrfl flPreibisch et al.|l999| |Schertl et aL|2003||Kraus et al. 
|2007 l, a first-order glance at the Trapezium shows the most 
massive stars in the cluster arranged in a central, compact 
configuration. Observations of a possible proto-trapezium, 
W3 IRS 5 (Megeath et al. 20051, seem to show a massive 



old enough to have segregated to the degree seen (Bonnell 
|fc Daviesj 19981 ). 



Recent work (IMcMillan et al. 20071 has suggested a 



subcluster still in the embedded phase, which could be sug- 
gestive of either formation as a compact cluster or gas-driven 
migration during formation. However, if a compact system 
of massive stars were to form in situ, it would be dynam- 
ically unstable ( Pflamm- Altenburg & Kroupa 20061; if the 



stars in the Trapezium formed in such a fashion, perhaps an 
order of magnitude more stars were originally in the system 
in order to leave behind the 4 OB stars today. Alternatively, 
the stars may have migrated there by some combination of 
dynamical mass segregation or gas dynamical effects during 
formation. The problem with dynamical segregation alone 



scenario in which merging smaller-n subclusters can create 
a situation in which the final larger-n cluster appears to 
be dynamically much older than the age of the stars would 
suggest. The key concept is that dynamical mass segrega- 
tion occurs on a shorter timescale in the low-n subclusters, 
and that the segregation is maintained in the merger pro- 
cess, so that effects of the shorter dynamical timescale of the 
subclusters are applied to the structure of the final merger 
products. This idea is supported by phase-space arguments 



by Fellhauer et al. (20091. This sort of later dynamical mix- 



ing is perhaps supported by observations suggesting that 
young massive stars are frequently found in close proximity 
to other stars of varying ages ( Brogan et al.|2008 i. 
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1 Though a complete theory for the formation of a Trapezium- 
like system needs to include this key aspect. 



The McMillan et al. results point the direction toward 
an attractive solution to the rapid concentration of massive 
stars as seen in the Trapezium and W3 IRS 5. However, their 
simulations employed large numbers of stars, with n = 10 4 - 
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4 x io 4 . it remains untestecj^] how applicable these results 
are to less prolific but more representative star clusters with 
n — 300-1000, which may be more representative of a 'typ- 
ical' star forming region ( Lada & Lada 20031 IPorras et al 



100 



2003). In these smaller clusters with grainier potentials and 
crossing times that are smaller than their larger siblings, it 
is unclear if the factor of ~ 2 decrease in the dynamical 
time gained by halving the number of stars will come into 
play. In this paper we investigate the merger of low-n sub- 
clusters to determine if this process is likely to play a role 
in accelerating the concentration of massive stars in young 
local clusters. 



2 MEASURING SEGREGATION 

There are several ways to measure the level of segregation 
in a cluster, including the evolution of the Lagrangian radii, 
the half-mass radii of different mass components, or their cu- 
mulative number distributions. A new method was recently 
developed by Allison et al. (2009a I, utilizing the minimum 



spanning tree (MST) of the stars of interest relative to that 
of a random subset of the cluster stars. An attractive feature 
of this method is that it doesn't rely on defining the center 
of a cluster. The method we use here is directly inspired by 
the MST method. 

Rather than use the MST, we use the two-dimensional 
convex hull of the stars' positions projected onto the plane. 
The convex hull encloses the minimal area that contains all 
straight lines joining any two points; in two dimensions it is 
conceptually found by stretching an elastic band around the 
points of interest and letting it relax, and it is unique for a 
given set of points. The area of the convex hull of N stars, 
which we denote H C (N), divided by N yields an effective 
surface density for the stars, and can be used like the MST 
to measure the compactness of a subset of stars in a cluster. 

Our scheme proceeds much as the MST method, in the 
following manner for the N most massive stars in the cluster: 
first, compute H C (N) for many subsets of N stars drawn 
randomly from the cluster. Computing the convex hull is 
not intensiv^] 10 4 subsets of 100 stars can be computed in 
~ 15 seconds on an aging Apple MacBook with a 2 GHz 
Intel Core 2 Duo and a loud cooling fan. After computing 
1000 or more values of H C (N), the median and ±ler, ±2<r, 
and ±3<r values of the distribution can be determined. These 
define the 'normal' range of values for H C (N) in the cluster 
for a random sampling of N stars. Second, compute H C (N) 
for the N most massive stars. Finally, repeat these steps 
for many values of N up to the highest interesting value. 
The compactness of the N most massive stars can then be 
compared to the expected value for that N, to see if they 
are more closely clustered than a random placing would be 
likely to yield. 

To illustrate this scheme, we use it to compute the segre- 
gation levels for three test scenarios: an unsegregated Plum- 



Although after submission of this paper, | Allison et at] ( |2009b| 
published a similar study to this one, which we address in the 

discussion. 

3 There are several algorithms for the calculation (e.g. |Press et al.| 
|2007| >, which in two dimensions scales as nlog(n) for n points. This 
work uses the built-in function implemented in Mathematica. 
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Figure 1. Examples of mass segregation as detected by the 
scheme used in this work. In each plot the black line is the median 
value of H C (N) for the given N, and the progressively lighter gray 
regions show the ±1<t, ±2(7, and ±3cr regions. The coloured lines 
show 10 realizations of a random setup for the Plummer and Subr 
clusters, and for the ONC case the data from [Hillenbrand ( 1997| ). 
The Plummer spheres and Subr clusters all have 1000 stars and 
half-mass radii ~ 0.5 pc. 
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Figure 2. Example of two merging clusters, each with n = 500. 
Red dots show the positions of the eight most massive stars in 
the experiment. 
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mer sphere; a segregated cluster using the method of Subr 
et al. (20081, using their value S = 0.5 corresponding to 



a dynamically mature cluster; and following Allison ct al 



( 2009a I , the Orion Nebula Cluster using the data from 



Hil- 



lenbrand (19971. The Plummer sphere and segregated clus- 
ters all have 1000 stars, and characteristic radii ~ 0.5 pc. 
These examples are shown in figure [l] In each panel, the 
black line is the median value of H C (N) for 3000 random 
samples of N stars, and the progressively lighter gray re- 
gions show the ±1<7, ±2cr, and ±3<r values of this distribu- 
tion. These regions are to be compared to H C (N) of the N 
most massive stars, which we plot as coloured lines. 

Looking at the Plummer sphere in the top panel, note 
that most of the coloured lines, showing the values of 10 in- 
dividual random realizations of the cluster setup, lie within 
the ±2<r shaded region of typical surface densities for each 
N. This is then the indication of an unsegregated cluster. 
Note that the median values and shaded regions are taken 
from a single realization; for these clusters and for the seg- 
regated cluster, the spatial distributions of each cluster are 
very similar and the median and shaded regions are almost 
indistinguishable between realizations. In the middle panel 
we show a heavily segregated cluster, and we see that for 
all cases H C (N) of the N most massive stars lies above the 
+2cr value up to N = 45, and for seven out of the ten this 
extends out past N = 100. This means that for most of 
the realizations of this setup, at least the 100 most massive 
stars are distributed in a fashion that is more compact than 
a random sampling of 100 stars would yield. 

Turning to the lower panel, we show the method as 



applied to the data from Hillenbrand ( 1997 1 for the ONC 



We see that for the lowest value of N that we consider, 
N = 4, the surface density of the most massive stars is quite 
high. For N < 10, H C (N) lies around or above the +3cr 
value of the cluster distribution, before falling back at or 
below the cluster median for higher values of N. This agrees 
quite well with the result | Allison et al.| |2009a) found for 
the same data, which is unsurprising given that this method 
is fundamentally indebted to theirs. Note also the higher 
values of the surface density for the ONC compared to the 
randomly generated clusters. This illustrates the retention 
of physical meaning this method offers, showing both the 
absolute surface density and the relative amount of mass 
segregation for the most massive stars in a cluster. 



3 NUMERICAL EXPERIMENTS 

We performed experiments to compare the amount of mass 
segregation obtained with a single Plummer sphere to that 
resulting from the merger of two smaller Plummer spheres. 
The systems were integrated using nbody6 ( |Aarseth|2003[ ). 
As these are merely illustrative experiments, we used a sim- 
ple one-component mass function with masses in the range 
0.1 - 10.0 MQand a power-law index a — —2.3. We ran cases 
with two final cluster sizes: larger clusters where the single 
cluster had n = 1000 and scale length a — 0.5 pc, and the 
merging clusters had n = 500 and a — 0.25 pc; and smaller 
clusters where the single sphere had n — 300 and a — 0.2 pc, 
and the merging clusters had n = 150 and a = 0.1 pc. These 
cluster sizes are a factor of a ~2-4 times more compact then 



et~aT1[2006t . If one considers the two merging spheres as a 
single cluster at a pre-merger time, however, they would be 
larger than that sample. This is also ignoring the (unmod- 
eled here) effects of gas dispersal and non-sphericity. We 
chose the upper mass limit of 10 Mq as we are interested 
primarily in low-to-moderate mass star formation. Addition- 
ally, our clusters are of relatively low mass (~ 100-325 Mq), 
where one would arguably not expect to find a star of much 
greater mass than 10 Mq ( Weidner &: Kroupa|2006 1. 

The merging clusters were created by setting up one 
Plummer sphere and setting its virial ratio, then duplicat- 
ing the cluster. The copies were initially at rest with respect 
to each other, placed so that they merged at an 'early' time, 
2.25 x 10 5 yr, and a 'late' time, 7.5 x 10 5 yr. Segregation 
was then measured at 3 x 10 5 yr and 10 6 yr for both the 
merger product and the single Plummer sphere. An exam- 
ple of this is shown in figure [2] As there is some evidence 
that star formation may result in sub-virial stellar veloci- 
ties ( e.g. |Walsh et aLl|2004| |Peretto et aT1|2006| |Kirk et al 



the clusters in Lada & Lada (2003 1 (see figure 2 of Adams 



20071, each setup was run from virialized and cold initial 
conditions, with virial parameters Q V i r = 0.5 and 0.1. Each 
set of parameters was run with 10 random realizations of 
the initial conditions. As a simple way to exclude escaping 
stars from the mass segregation calculation, we restrict our 
analysis of the results to stars that remain within 1.5 pc of 
the center of mass. 

We have anchored our experiments in physical time 
from the commencement of the integrations, rather than in 
terms of the initial crossing time, which only has observa- 
tional significance to the single Plummer spheres with their 
relatively static structure. We therefore also determine for 
each set of parameters an apparent dynamical age, which is 
what a present-day observer might measure (assuming rea- 
sonably accurate stellar ages can be determined). We cal- 
culate this by determining the crossing time at the end of 
the experiment and dividing this by the age of the sys- 
tem in years, i.e. T dy „ = t V r/t cro ss, final- We define the 
crossing time as t cro3S = 2r; l /u rms , with ru the half-mass 
radius and v rms the root mean squared velocity. We also 
show the age of the system in terms of its final relaxation 
time, Trelax — tyr j trelax , final, using the standard calculation 
for the relaxation time t re i ax = N / (8\nN)t cross . The mean 
value for each of these quantities for all sets of parameters 
is shown in table 1. 

The results of these experiments are shown in figures 
[3] [6] where we plot the segregation up to the 60 most mas- 
sive stars for each series of runs. Note that in these plots the 
median and its spread are calculated for a single realization, 
and all ten individual realizations are compared to this one. 
While there are small differences in the median distribution 
between runs, for clarity we show only one, and this does not 
affect our conclusions. In table 1, we summarize these plots 
and show the fraction of runs (out of 10) Ni that display seg- 
regation of the i most massive stars, for i — 5, 10, 20, 40. For 
these purposes we define 'displaying segregation' as having 
a value of H C (N) above the +2a value for the typical cluster 
value. We also plot these fractions for the n = 1000 cluster 
in figure [7| Examining these figures and the table, it is clear 
that within a pair of runs with the same age, n, and Q v i r the 
cluster resulting from a merger displays more segregation, 
with the difference between solitary Plummer spheres and 
the merger products being more marked among the higher 
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Figure 3. Mass segregation from our runs with 1000 stars total, with virialized initial conditions. As in figure [I] the coloured lines are 
10 random realizations of the initial conditions, and the gray shaded regions show the ±1, ±2, and ±3 sigma values of H C (N) about the 
median black line. Time in years and the apparent dynamical age are labeled. 
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Figure 4. As figure [3] for subvirial initial velocities. 
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Figure 5. Mass segregation from our runs with 300 stars total, with virialized initial conditions. As in figure [l] the coloured lines are 
fO random realizations of the initial conditions, and the gray shaded regions show the ±1,±2, and ±3 sigma values of H C (N) about the 
median black line. Time in years and the apparent dynamical age are labeled. 
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Figure 6. As figure[5] for subvirial initial velocities. 
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Figure 7. The fraction of runs A/j showing segregation for the i most massive stars, for the n = 1000 clusters. 



n clusters. As we have chosen to fix the integration time to 
physical times rather than initial crossing times, this is un- 
surprising, since the lower- n clusters are dynamically more 
mature. The difference between the Plummer spheres and 
the merger products is less pronounced with sub-virial ini- 
tial conditions, and interestingly the n — 1000 and n = 300 
sub-virial Plummer spheres have identical apparent dynam- 
ical ages at the end of the simulations. 

It is also apparent that none of the initial conditions 
we examine here lead to a cluster displaying the full segre- 
gation of the clusters created using the technique of Subr 



et al. (20081, which corresponds to the end-state of dynam- 
ical mass segregation. The segregation in these simulations 
dominantly affects the most massive stars, a situation that 
looks more like the ONC data from |Hillenbrand| ( |1997| ) in 
figure [l These results would then seem to confirm the find- 
ings of McMillan et al. (20071 in the lower- n regime that 



may be more applicable to local star formation. As the early 
evolution of young clusters with imposed primordial segrega- 
tion across all masses appears to be inconsistent with young 
observed systems ( jMoeckel &: Bonn ell 2009 ) , this scenario is 
attractive since it seems to affect only the most massive stars 
in a cluster over these timescales and population levels. A fi- 
nal point is that the single Plummer sphere models with cold 
initial conditions exhibit more segregation of the most mas- 
sive objects than their virialised counterparts; while not at 
the same level as the merged clusters, this is another possible 
route to the formation of dense Trapeziumesque systems. 

Examining the pairs of runs of the same physical age in 
terms of their apparent dynamical age, the general trend is 
that the merger product appears to be 1-4 dynamical times 
more advanced than the solitary Plummer sphere. The ex- 
ceptions to this are the late-merging sub-virial case with 
n = 1000 and the late-merging virial case with n = 300, 
both of which appear to be dynamically younger than their 
solitary counterpart, though the difference is very slight. 
The sub-virial initial conditions for the n — 1000 Plummer 
spheres actually appear to be dynamically older than their 
virial, merged counterparts, showing that both merging and 
sub-virial initial conditions can make a cluster appear older 
than its true age The fraction of runs showing segregation is 
not directly related to the apparent dynamical age, however; 
the snapshot afforded by an observation may not reveal the 
history that led to that point. In the sub-virial Plummer 
sphere runs, for instance, the evolution of Td yn rises sharply 



from to ~ 2.5 and plateaus for several hundred thousand 
years, as the cluster collapses, overshoots virial equilibrium, 
and expands. 

Considering the apparent relaxation times at the end of 
the simulations, one would not expect complete segregation 
as a result of two-body relaxation, as none of the systems are 
older than trelax- Since the segregation timescale for stars 
of mass m decreases as < m > /m, however ( Spitzer|1969 1, 
one would expect to see some segregation in the massive 
stars based on the apparent relaxation times. The fact that, 
for instance, the virial Plummer sphere with n — 1000 and 
its late-merging counterpart have very similar ages relative 
to the relaxation time, yet very different mass segregation 
characteristics, highlights the limited applicability of the re- 
laxation time in young systems of unknown history. In more 
complicated star formation scenarios (i.e. non-spherical, or 
out of virial equilibrium), observational estimates of the 
crossing time or relaxation timescale should be treated with 
caution, as also emphasized by Bastian et al. ( 2008[ ) and 



Allison et al. ( 2009b I 



4 DISCUSSION 

A cautionary note should be sounded about the interpreta- 
tion of these results and similar early-cluster stellar dynam- 
ical experiments. There have been several attempts in the 
literature to apply n-body techniques to very young clusters, 
when a proper treatment should necessarily include the gas 
reservoirs out of which the stars are born. To the extent 
that the remant gas is treated, it is usually as a spherical 
or spheroidal background potential on the cluster scale (e.g. 
|Adams et al.|2007| [Proszkow et al.|2009[ ), or as an attempt 
to interpolate simulated hydrodynamic results on the level 
of individual stars (Moeckel & Bally 2007). As it becomes 
more clear from observations and simulations that cluster 
formation may not be well-described as a spherical process 



(e.g. |Myers et al.|199l] |Bonnell et al.|2008| |Proszkow et al 



2009), work from the n-body end has pushed toward clus- 
tered and asymmetrical initial conditions as in |McMillan| 



et al. (20071. It is not immediately clear that divorcing the 
stellar dynamics from the gas dynamics and treating effects 
like mass segregation as purely dynamical processes should 
be succesful; the inclusion of gas dynamics as in the large- 
scale simulations as in, for example, |Bonnell et al.| ( |2008[ ) or 
Bate (2009) may prove to be obligatory. 
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The background potential of any remaining gas can be 
included in an n-body simulation with relative ease. Two 
immediately relevant dissipative effects that cannot be mod- 
eled, but are present in hydrodynamic simulations, are the 
gas drag acting on the stars as they orbit in the cluster 
potential, and the mass loading that occurs in a competi- 
tive accretion scenario. Both of these effects scale up with 
the mass of the star, and both should act to accelerate mass 
segregation. Interestingly, |Bate] ( |2009| reported minimal seg- 
regatiorj^Jin his recent simulation that included these effects 
as well as merging sub-clusters, indicating that (as is usu- 
ally the case) matters are more complicated than simple 
prescriptions might allow for. This suggests that while ef- 
fects like deceleration due to Bondi-Hoyle accretion could 
perhaps be added to an n-body code, the utility of such 
approximations may be limited. 

A comparison between pure n-body merging experi- 
ments can be made to the recent hydrodynamic/n-body 
simulation presented in Bate (2009). In this simulation, the 
cluster at 2.85 x 10 J yr is the end result of the merger of 
five sub-clusters consisting of ~ 1250 stars. These mergers 
took place in the presence of the cluster gas. Bate notes that 
mass segregation is not present in the end cluster as mea- 
sured by the cumulative radial distribution, or the median 
mass, upper quartile mass, and maximum mass except for 
in the inner 1000 AU. Applying the segregation detection 
method used in this paper to the positions of the stars in 
Bate's simulation, we find evidence of segregation for the 
most massive stars, shown in figure [8] While there is clearly 
not general segregation as in the Subr clusters, the ~ 10 
most massive stars display a level of segregation similar in 
character to the ONC data, and would not stand out as out- 
of-place in amongst the simulated data in figures [3] and [4] 
While the issue of the equivalence between gas-free and em- 
bedded stellar dynamics should be explored further, in this 
case the level of segregation seen appears to be similar in 
character. 

The mass segregation signal resulting from the merger 
of sub-clumps of an age and stellar content suitable to 
nearby young clusters appears to affect predominantly the 
most massive stars in the cluster, and does not present it- 
self as a general trend that might appear in a cumulative 
distribution analysis or similar techniques. Full-scale hydro- 
dynamic simulations ( |Bate|2009[ ) that feature merging sub- 
clusters display a similar signal, despite the mergers tak- 
ing place in the presence of gas. While the more general 
type of segregation may have been looked for as the 'goal' 
of segregation-producing phenomena in the past, the type 
that affects only the most massive stars may be more reflec- 
tive of Nature (as seen in the ONC), and thus sub-cluster 
merging and perhaps sub-virial initial conditions seem to be 
a promising scenario. 



1000 stars placed via a fractal distribution. The rapid seg- 



regation seen in the Allison et al. ( 2009b I experiments is 



Allison et al. ( 2009b \ published the results of an inde- 



pendent investigation that overlaps significantly with this 
work. Namely, they studied the dynamical effects of sub- 
virial, clumpy collapse on the mass segregation in a cluster. 
Their initial cluster is a more naturalistic setup versus the 
stylized two-Plummer-sphere clusters we consider here, with 



4 Though minimal, measurable in the innermost regions; see be- 
low. 



mostly due to two-body relaxation in the dense core that 
develops as the sub-clusters merge and the cluster collapses 
at early times. The extent to which pre-merger segregation 
in the sub-clusters contributes to the final state (McMillan 
et al. 2007; Fellhauer et al. 2009; this work, see figure [5| 
versus two-body relaxation in the dense core resulting from 



clustered, sub-virial collapse (Allison et al. 2009b) is likely 
dominated by the details of the initial setup. 

Finally, we caution against an overly-broad definition of 
Trapezium-like or Trapeziumesque. Figure [9] shows the lo- 
cations and clustering dendrogram of the four OB systems 
in the Trapezium and their entourage ( Preibisch et al.|1999 



Schertl et al.|2003|[Kraus et al.|2007| ). Most work, including 
this one, focuses on the four main branches of the dendro- 
gram, i.e. attempting to create a clustering of four or more 
massive stars. While this is clearly the first step toward a 
Trapezium, there are at least two more levels of clustering, 
in the mini-Trapezium of O x B and the wider binary of the 
O A hierarchy, and then the closer visual and spectroscopic 
binaries of all four main components. Turning to the younger 
W3 IRS 5, with the further complication of surrounding na- 
tal gas, there is evidence for multiplicity at sub-1000 AU 
scales (Megeath et al. 20051. Bearing in mind the limited 



number of examples, multiplicity of the massive sources in 
Trapezium-like systems seems ubiquitous, and could be fun- 
damental rather than a detail to be worried about later. If 
dynamical effects like those considered here and in |Alliso"n| 
et al. ( 2009b I are to explain the clustering of massive stars, 



these processes must either leave intact the presumably pri- 
mordial multiplicity of the massive stars, including higher 
order subsystems like 1 B, or produce the multiplicity as 
part of the segregation; neither possibility has yet been in- 
vestigated, and consistency with this aspect of the Trapez- 
ium may help discriminate between dynamical segregation 
versus birth location as the cause of massive star clustering. 
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Table 1. Cluster characteristics at the end of the integration. All mean quantities are measured at the end of the simulation. 
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Figure 8. Mass segregation from the end state of the cluster 
presented in |Bate| ( [2009] ) . 
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Figure 9. The locations and hierarchy of the Trapezium stars. 
The length scales assume a distance of 400 pc. 
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